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Abstract 



It has been empirically observed that the dependence of J /ip and ifj' production 
on nuclear mass number A is very similar. This has been postulated to be due 
10 '. to the predominance of color octet pre-resonant states in charmonium produc- 

er I tion and absorption. Two new experiments, NA60 at CERN and HERA-B at 

r~^ ■ DESY, will measure the Xc A dependence for the first time. These measure- 

^^ , ments should shed new light on the charmonium production and absorption 

^^ ' mechanisms. 
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1 Introduction 

The description of particle production in proton-nucleus interactions as a simple scal- 
ing of the proton-proton production cross section with nuclear mass number A, 

(^pA = CTpNA."' , (1) 

has been used to describe many processes. When the production of the desired 
final state particle is calculable in perturbative QCD, the factorization theorem |]I| 
suggests that the exponent a should be unity. Drell-Yan production, integrated over 
all kinematic variables, agrees with a = 1 to rather high precision although some 
deviation from unity appears at large values of Feynman x, xp = p\\/pma.x- A less 
than linear A dependence has been observed for J/ip, '^', TT, and T' + T" production 
with 0.9 < a < 1 near xp ~ 0. 

By now, the A dependence of J/tp production at Xi? > is known to rather high 
precision at several different energies |^, ^, |^, |^, 0, |^ . While the ip' A dependence is 
not as accurately known, its statistics were sufficient for the E866 collaboration to de- 
termine that the ip' a is smaller than the J/ip a for xp < 0.2 P). The known A depen- 
dence of J/ip production has been used to determine the strength of the "anomalous" 
J/ip suppression in Pb+Pb interactions at the CERN SPS [^. However, an important 
assumption in this interpretation is that all charmonium states interact with the nu- 
cleus while in "pre-resonant" \{cc)sg) states [0. Since a significant fraction, ~ 40% 



0, 



of the observed J/ip^s come from Xc decays [0, |T2|, |T3|, |T^, [T^, 0, |T^ |T8|, |T^, |2^, ^ , a 
measurement of the Xc A dependence is crucial for the understanding of J/ip suppres- 
sion in nucleus-nucleus collisions because in a quark-gluon plasma, J/ip suppression 



is expected to occur in steps, the first of which is the dissociation of the Xc [H] • So 



far, only one experiment [^ has presented differential distributions of Xc production, 
albeit with such a small sample that it is not possible to tell if the shapes of the Xc and 
J/ijj xp distributions are the same or not. No measurement of the Xc A dependence 
has yet been made. 

Fortunately, this situation seems about to change. The Xc A dependence will 
be measured for the first time in two fixed-target experiments at different energies. 
The NA60 collaboration, a follow-up to the NA50 collaboration at CERN, has been 
approved for pA measurements at 450 GeV and is planning to also take data in 



nucleus- nucleus interactions at 158 GeV p3|, |2^. Their muon spectrometer will sit 



at < i/cm < 1 at both energies, giving forward xp coverage only. The HERA-B 
collaboration at DESY has placed target wires around the halo of the proton beam 
at HERA. In their first run, they have demonstrated that they can detect the X' 



c 



25| , p6| . In their next run, they will measure the J/^, ^', and Xc A dependence over 
—0.5 < Xp < 0.3. This will thus be the first charmonium experiment with coverage 
significantly below xp ~ —0.1. 

If the A dependence of Xc production is the same as that of the J/ip, then the 
picture of a pre-resonant color octet state passing through the target |T0| will be 
validated. Then charmonium production and absorption at fixed-target energies can 
be essentially described within the color evaporation model [E^]. However, the non- 



relativistic QCD model |28, 29 1 predicts that Xc production should be predominantly 



color singlet while direct J/ip and ip' production is via color octet states. If this pic- 
ture is correct, the A dependence of Xc production could be quantitatively different 
than that of the J/ip and ip' . 

In this paper we focus only on charmonium production and its subsequent ab- 
sorption by nucleons, the "normal absorption" identified by NA50 ||^. While this is 
insufficient to describe a{xF) over the full range of xp-, the NA60 and HERA-B mea- 
surements will be in a region where the xp dependence of a has so far either not been 
determined or has not been strong ||^. At larger negative xp-, the A dependence may 
be different than expected from pre- resonant absorption |TD|. Other nuclear effects 
such as shadowing, energy loss, and intrinsic heavy quarks depend only on either the 
projectile or target momentum fractions and not on the identity of the final charmo- 



nium state ||3^ and thus should affect J/ip and Xc production identically. We first 
discuss charmonium production by color evaporation and nonrelativistic QCD and 
then describe how nuclear absorption of color octet and color singlet states might be 
disentangled. 

2 Charmonium Production: Color Evaporation vs. 
NRQCD 

Two models have been used to describe quarkonium hadroproduction: the color evap- 
oration model (CEM) and the nonrelativistic QCD model (NRQCD). Since both have 
been described in detail elsewhere, we only discuss the specifics that are germane to 
our calculation. 

In the CEM, charmonium production is a subset of the cc pairs produced below the 
DD threshold. The hadronization of charmonium state C from these sub-threshold cc 
pairs is accomplished through the emission of one or more soft gluons. It is assumed 
that the hadronization does not affect the kinematics of the parent cc pair so that 
only a single universal factor, Fc, is necessary for each state. The factor Fc depends 
on the charm quark mass, mc-, the scale /i of the strong coupling constant a^, and 



the parton densities. We use the MRST LO parton distributions |^, ^ for CEM 
production. The factor Fc must be constant for the model to have any predictive 
power. The differential and integrated quarkonium production rates should thus be 
proportional to each other and independent of projectile, target, and energy. The 
relative charmonium rates so far seem to bear this out since Yl,j Xcj/{J/'4') ~ 0-4 and 



ijj' /{J /ip) ~ 0.14 ||T8|, [T^, ^ |33|, |3^ over a wide range of targets and energies, see also 
Ref. [gg. 

The LO cross section of state C, ac-, from projectile p and target A is 

- — = 2F(j K / mdm / dxidx26{xiX2S — m^) 5{xf ~ Xi + X2) 
dxp J2mc Jo 

X l^fP{xi,m^)ff{x2,m^)agg{m^) 

+ Y: [f^ix^,m')f^^ix2,m') + f^ix^,m')f^ix2,m')]a,,{m')] . (2) 



= U,d,! 



The partonic cross sections (Tgg and aqq can be found in Ref. |3^, |3^. Production by 
quark- gluon scattering enters only at NLO. 

A K factor was included in Eq. (^ since our calculation is at leading order and Fq 
was determined at next-to-leading order, as indicated. At NLO, the charmonium cross 



section was calculated using the QQ production code of Ref. |^ with a cut on the pair 



mass as in Eq. (H) |27|. The pt dependence and the normalization of the charmonium 



cross section from the Tevatron collider agrees with these calculations [^]. Since, at 
fixed energy, the K factor for cc production is independent of the kinematic variables 
, our calculation is at leading order. Therefore, we multiply the LO cross section 



by K to obtain the magnitude of the NLO cross section and then also by Fq^'~' to 
fix the hadronization of the subthreshold cc pairs to charmonium. Note however that 
since we study ratios of cross sections, only the relative normalization is important and 
because no nuclear effects on the parton densities are included, the GEM production 
information generally cancels. 

Since the CEM depends on the universality of charmonium hadronization through 
soft gluon emission, a check of this assumption for Xc production, particularly as a 
function of xp, is critical. The Xc has previously been crucial for furthering the 
understanding of charmonium production. The color singlet model (CSM) [BO, ET 



described high px charmonium production as direct color singlet production with the 
appropriate quantum numbers. In the CSM, direct J/ijj and tjj' production required 
the emission of a hard gluon and should thus be rare on a perturbative timescale. 
However, x^'s could be directly produced as color singlets and thus high pt J/ip Pro- 
duction should be dominated by Xc decays. The measurement of Xc relative to direct 
J/ip production at the Tevatron collider ^^ showed that the CSM was incomplete. 

The non-relativistic QCD approach to quarkonium production was formulated 
p8|l was formulated as a way to go beyond the CSM. NRQCD describes quarkonium 
production as an expansion in powers off, the relative Q-Q velocity. Thus the angular 
momentum or color of the quarkonium state is not restricted to only the leading color 
singlet state but includes color octet production as well. 

The xp distribution of charmonium state C in NRQCD is 



dac 



dxp 



^^ f'dx,dx2S{xF-X,+X2)fr{x„fi')ff{x2,fi')C%^^^{0^) , (3) 



where the partonic cross section is the product of perturbative expansion coefficients, 
C*^r ,, and nonperturbative parameters describing the hadronization, (O^). We use 
the parameters determined by Beneke and Rothstein for fixed-target hadropro duct ion 
using the CTEQ 3L parton densities |^ with iric = L5 GeV and /i = 2mc ||29|. Since 



the parameters {O*^) are fit to the LO calculation with a LO set of parton densities, 
no further K factor is required. 

Direct J /xp production has only contributions from gg fusion and qq annihilation 



29|| , as in the CEM. The qq contribution is all octet while the gg component is a 
combination of octet and singlet production. The gg partonic cross sections for J /%p 
and %!)' production are 

d{gg ^ xlj) = Cll,AOt) = B^{x,,X2,s,m^)AM) + B,{x,,x^,s,ml){Ot{'S^)) (4) 
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Table 1: The percentage of charmonium production from color octets in NRQCD at 
each energy we consider. 



where only the octet, Af = (O^(^5o)) + (7/m2)(0|'(3Po)), and singlet, (Of (^^i)), 
matrix elements differ between J /ip and ip' production. The functions Bi and B^ are 
proportional to a1 and a^ respectively. The octet parameters are quite different for 
the two states: As{J/ip) ^ 5.8A8(?/''). The smaller Agi^ip') could be due to the larger 
mass and thus the increased "hardness" of the emitted gluon for the ip'- 
On the other hand, a color singlet Xc can be formed from two gluons 



40, 41 



so 



that Xc production is predominantly color singlet. In addition, the Xci has a singlet 
contribution from gq scattering at 0{a^) [^. Only the qq channel contributes to color 
octet Xc production. Thus the largest singlet contribution to total J/ip production 
is from Xcj decays. Of these Xcj decays, the most important is the xa which has a 
27% branching ratio to J/ip. The Xc2 also has a relatively large branching ratio to 
J/ip, 14%. Although the Xco production cross section is as large as those of the other 
Xc states, its small branching ratio, < 1%, results in a negligible Xco contribution 
to J/i/j production. The Xco is essentially invisible in hadroproduction experiments 
which reconstruct Xcj's from their radiative decays to J/ip. 

The total J/ip xp distribution then includes radiative decays of the three Xcj 
states and hadronic decays of the ■?/'', 



da 



J/^ 



dxi 



dxp 



+ E5(x 



J=0 



dxp 



J/i,X) 



da. 



jjj' 



dxi 



(5) 



In Fig. |1] we show an example of the relative singlet and octet contributions to total 
J/ip, direct J/ip, ip' and the sum of the three Xc contributions to J/ip production 
at 450 GeV, the SPS proton beam energy. Only the forward xp distributions are 
shown since the distributions are symmetric around xp = 0- No nuclear effects on 
the parton distribution functions are included. 

The percentage octet production of each charmonium state is given in Table |I|. 
The octet contribution decreases with energy for all charmonium states. Since color 
singlet Xcj production is through the gg and gq channels, the fraction of octet Xcj 
production is quite small, 10% or less. The decrease of octet production with energy is 
expected because the octet qq contribution becomes even smaller at higher energies. 
On the other hand, direct J/ip and ip' production is color octet dominated since 
both the gg and qq channels have octet contributions, see Eq. (^). The larger value 
of A^iJ/ip) relative to Ag (■?/'') increases the overall octet contribution from ~ 66% 
for the ip' to ~ 87% for the direct J/ip at ^/s = 29.1 GeV. However, when the 
XcJ radiative decays are included, the octet contribution to total J/ip production is 
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Figure 1: Charmonium xp distributions at 450 GeV. The total J/ip (a), direct J/ip 
(b), ip' (c) and summed Xcj contributions to the J/ip (d) cross sections are shown. 
The octet (dashed) and singlet (dot-dashed) contributions to the total (solid) are 
shown separately. 



nearer to that of the ip' , ~ 63%. These results are reflected in Fig. |l]. Note also 
that the xp distributions of the charmonium states are not exactly parallel to each 
other, as predicted by the CEM. Unfortunately the slopes of the xp distributions 
are quite similar and it is not until relatively large values of xp that the differences 
become more significant. However, there are other ways to distinguish the production 
mechanism since these two models of charmonium production lead to quite different 
predictions of the A dependence, as we will demonstrate in the next section. 

3 Absorption by Nucleons 

In Ref. [1^, absorption was described in terms of the singlet and octet components 
of the J /ij) wavefunction, 

\J/i)) = aol(cc)i) + ai\{cc)8g) + a2\icc)igg) + a2\{cc)8gg) H . (6) 

In the CSM [^, ^, only the first component is nonzero for direct J/ip production. 
The cc pairs then pass through nuclear matter in small color singlet states and reach 
their final state size outside the nucleus, at least when x^ > 0. If cc pairs are 
predominantly produced in color octet states, then it is the |(cc)85') state that interacts 
with nucleons. After the color octet cc is produced, it can neutralize its color by a 
nonperturbative interaction with a gluon. This octet state is fragile so that a gluon 
exchanged between it and a nucleon would separate the (cc)8 from the gluon, exposing 
its color and, since the octet is unbound, break it up [|1^. If the |(cc)85') state is free 
to evolve without interaction, such as in pp collisions, the additional gluon would 
be absorbed by the octet cc pair, hence 'evaporating' the color. The CEM does not 
then care about the relative coefficients in Eq. (P). As formulated in Ref. |^, the 



NRQCD model provides the leading coefficients in the expansion of the wavefunction 
in Eq. (|^) and hence encompasses both singlet and octet production and absorption. 
In this section, we will describe the absorption of color singlets, color octets, and the 
combination of the two for final-state J/ip, ip' and Xc production. Any differences in 
the A dependence of these states will be a consequence of this nucleon absorption. 
We calculate charmonium production in the CEM with pure octet and pure singlet 
absorption while NRQCD is used to determine the fraction of charmonium states 
production in color singlets and color octets. This then determines the rate of singlet 
and octet absorption in Eq. @. 

The effect of nuclear absorption alone on the J/ip production cross section in pA 
collisions may be expressed as |Q 

//•oo 
d% dzpA{b,z)S'''''{b,z) (7) 

J —oo 

where b is the impact parameter and z is the longitudinal production point. When 
the production and absorption can be factorized, as in the CEM, and no other A de- 
pendent effects are included, ap^ is independent of A and drops out of the calculation 
of a. The nuclear absorption survival probability, S'^^^, is 

5-^^(6, z) = exp {- ^" dz'pAib, z')a^Uz' - ^)} (8) 



The nucleoli absorption cross section, o"abs, depends on where the state is produced 
and how far it travels through nuclear matter. Nuclear charge density distributions 
from data are used for pA |0 • The effective A dependence is obtained from Eqs. (^ 



and (P) by integrating over z', z, and h. The full dependence on A can be related 
to a(xi?) in Eq. (|1|) but a is only constant if (Tabs is constant and independent of 
the production mechanism [^, Q. The observed J/ip yield includes an ^ 30% 



contribution from Xcj decays |5^ and an ^ 12% contribution from Tp' decays 
Then the total J/tp survival probability is 

Sf/lib, z) = 0.585}^;, dir(&, ^) + 0.35j^(6, z) + 0.12S;Y{b, z) . (9) 

The ip' and Xc states are only produced directly since other, more massive, charmo- 
nium resonances lie above the DD threshold and decay to DD pairs. 

We will present calculations for the total and direct J /ip, ip', and Xcj -^ J/'P A 
dependence. We include the Xcj branching ratios to J/ip because even though the 
Xco cross section is large, the small branching ratio gives it a negligible contribution 
to the final-state J/ijj yield. Our results will be calculated at 158, 450, and 920 GeV, 
corresponding to the NA60 and HERA-B energies respectively. We calculate a{xF) 
for two targets in each experiment: Be and Pb for NA60; C and W for HERA-B. 

3.1 CEM: color singlet absorption 

We first discuss pure color singlet absorption. In this case, (Tabs depends on the size 
of the cc pair as it traverses the nucleus. This was first described in terms of color 
transparency [Q . The cc pairs are initially produced with a size on the order of its 



production time, ri^a ~ Tjnit oc m^^. This initial size is ignored in the calculation. 
The charmonium formation time obtained from potential models [^ is re ~ 1 — 2 



fm, considerably longer. The absorption cross section of these small color singlet pairs 
grows as a function of proper time until tc when it saturates at its asymptotic value 

1471, ra, 



'^CN 



T 



aabs(.'-.)= ""H^J '^"<"^ . (10) 

[ (T^^ otherwise 

The proper time r is related to the path length through nuclear matter by r = 
{z' — z)/'yv where the 7 factor introduces xp and energy dependencies to (Tabs- At low 
energies and negative xp, the cc pair may hadronize inside a large nucleus. 

Figure illustrates the energy dependence of color singlet absorption in pA in- 
teractions. We take ctji^j^ = 2.5 mb |^. Assuming that the asymptotic absorption 
cross sections scale in proportion to the squares of the charmonium radii ED 



we 



have (T^/jv ~ 3.7(t}/^^ and ct^^^at ~ 2.4crj/ .^y. Thus each contribution to Eq. (H) 
has a different A dependence. The charmonium formation times are: tj/^ = 0.9 fm, 
r^/ = 1.5 fm, and r^^ = 2.0 fm [^. The results at xj? < reflect the differences in 



formation times as well as the change in the 7 factor due to their masses. The shorter 
formation time of the J/ip allows it to reach its asymptotic size at large negative xp- 

8 



At the lowest energy, 158 GeV, the charmonium states have a small chance of be- 
ing formed inside the target at Xi? > since a 7^ 1 at Xi? ~ 0.2 although the deviation 
from unity is small. For higher energies, the charmonium states are produced outside 
the nucleus for Xi;' > so that a ~ 1. No observable differences appear between 
the charmonium states at positive xp- Indeed, this A dependence is in contradic- 
tion with all available data at Xi;' ~ unless other nuclear effects are included |30|. 
Therefore this picture of absorption is primarily useful for the interpretation of our 
calculations of pure octet absorption in the GEM and the combination of singlet and 
octet production and absorption in NRQCD. 
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Figure 2: The A dependence for color singlet absorption is shown. The results are 
calculated at 158 GeV (a), 450 GeV (b), and 920 GeV (c). The total Jj^ (sohd), 
the direct J j'^ (dashed), the ■?/'' (dot-dashed) and the Xc (dotted) A dependencies are 
given. 

The direct J j^^ A dependence (dashed curve) is weakest because its asymptotic 
cross section is smallest. The ■?/'' A dependence (dot-dashed curve) is strongest because 



its final-state size and corresponding (Jabs is largest. In the calculation of Ref. [22 



the Xc radius is somewhat smaller than that of the ?/'' so that crLjy > ^% n- The x^ 



'XcN- 



formation time is the longest of the charmonium states and thus most likely to be 



produced outside the target. Therefore the Xc « is actually slightly larger than that 
of the direct J /ip aX xp ^ ^ due to the longer Xc formation time (dotted curve). The 
Xc contribution to the A dependence of the total J /ip yield (solid curve) decreases 
the total J /ip a at large negative xp, more like the Xc while when Xi? — »• 0, the Xc ot 
is near unity and the total and direct J/ip A dependencies are the same. 

3.2 CEM: color octet absorption 

On the other hand, if the cc pairs are produced only in color octet states, they should 
hadronize after rg ~ 0.25 fm in the cc rest frame ^9|. At large xp in the lab frame. 



hadronization then occurs after the cc has passed through the target as an octet. 
These fast cc pairs thus remain color octets until after they have left the nucleus. 
However, at negative xp it is possible for the octet states to neutralize their color 
inside the nucleus and interact as color singlets during the remainder of their path 
through the target [^]. This effect has typically been neglected when studying the 
A dependence of quarkonium production because the effect remains small in the xp 
regions so far covered, xp > —0.1 fS^, Q. (See however Ref. |^.) While travehng 



through the nucleus as a pre-resonant \{cc)8g) state, the eventual identity of the final- 
state resonance is undetermined and all quarkonium states are absorbed with the 
same cross section, a^^^. This physical picture agrees rather well with the empirical 
evidence that the J/ip and ip' A dependencies are similar over the measured xp range 
P, If. We choose cr°|-,g = 3 mb to agree with a ~ 0.95 for the J/ip measured by the 
E866 collaboration at Xi? = |p when no other nuclear effects are considered. Note 
that this value is somewhat smaller than typically assumed for the color octet cross 



section [T^] due to the relatively large measured a. 

We account for color neutralization of the octet in the nucleus in the following 
way: The path length of the Kccjsg) through the nucleus is calculated in the nuclear 
rest frame. If it exceeds the maximum path length through the nucleus from the 
\{cc)sg) production point, cr°bs = 3 mb for all charmonium states. This is the case 
for Xp > with all three energies. However, if color neutralization occurs before the 
state escapes the target, the resulting color singlet is absorbed according to Eq. (0). 

The A dependence of color octet absorption is shown in Fig. ^. Note that at 158 
GeV, color neutralization is achieved for xp < —0.2. At higher energies, neutraliza- 
tion occurs in the target at larger negative xp. It is important to remember that 
just because the octet color has been neutralized, the asymptotic cross section is not 
necessarily reached inside the target. With a formation time of less than 1 fm, only 
the J/ip is likely to be fully formed in a large nucleus, as observed in the 'saturation' 
of the A dependence at xp < —0.5. On the other hand, although the ip' and Xc niay 
become singlets inside the target, they do not reach their final-state size inside the 
target, even at xp -^ —1, due to the combination of their larger radii and longer 
formation times. The ip' and Xc A dependencies thus do not saturate, even at low 
energies. The slightly higher a at 920 GeV is due to the different target A ratios 
chosen for the NA60 and HERA-B calculations, Pb/Be and W/C respectively. 

We point out that the calculated a is lower in the pure color octet picture at 
Xp -^ —1 than in the color singlet absorption model even though the same asymptotic 
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Figure 3: The A dependence for color octet absorption is shown. The results are 
calculated at 158 GeV (a), 450 GeV (b), and 920 GeV (c). The total J/V' (sohd), 
the direct J j^^ (dashed), the V'' (dot-dashed) and the Xc (dotted) A dependencies are 
given. 

color singlet cross sections are used. This is because now the state starts out as a 
color octet with a finite probability to be absorbed before neutralizing its color. The 
probability tends to be larger in real nuclei where the path length is calculated in the 
integral over impact parameter rather than in the empirical analytic model where an 
average path length is used [E^ . The color octet is absorbed with its full cross section 



which is larger than the color singlet cross section at the point of absorption even 
though the asymptotic color singlet cross sections may be greater, e.g. a^,^ < a^/jy- 
The effective octet absorption cross section is larger because in the color octet state 
absorption has essentially no time delay. 

The greatest differences in the A dependencies of the states are at intermediate 
to large negative xp and would be most easily observable by NA60 at 158 GeV if 
their coverage extended so far. Note that, in this case also, the total J/?/^ and the 
Xc A dependencies would be quite similar in the target region while the ip' a would 
be lower, ~ 0.9 at 158 GeV and xp = —0.5 compared to ~ 0.93 for the total J/ip 



11 



and Xc- However, only HERA-B has the capabihty to measure the A dependence at 
negative xp and at the higher energy the differences appear at higher negative xp and 
are generally not as large. Thus any distinction will be rather difficult to determine 
and the observed A dependence is likely to be the same within the experimental 
uncertainties for all charmonium states. 

3.3 NRQCD: color singlet and color octet absorption 

Recall that in the preceding discussion, only one type of cc color state is assumed 
to be produced, either singlet or octet. Therefore the absorption factorizes from the 
production mechanism and the CEM cross section cancels in the calculation of a, as 
in Eq. (0). The result is then independent of all parameters in the production process 
such as TTic and the parton densities. However, according to Eq. (^, charmonium pro- 
duction is through a combination of octet and singlet states. In this case, production 
and absorption are intimately related and the NRQCD cross section determines the 



relative octet proportion for each state as a function of xp 0, |5^- The ratio of octet 
to singlet production is energy and xp dependent, as shown in Table ^ and Fig. |I]. 
Therefore Eq. (|^) does not hold since ap^ and S^^^ do not factorize for all xp. 

We now give the unfactorized xp distributions for each state in NRQCD. The xp 
dependence of direct charmonium production and absorption is straightforward: 
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(12) 



where ip = J/ip, ip' and T^ = J dzpAS^^^ for both singlet and octet absorption. The 
pp subscript is used to denote unmodified parton distributions in the target. The 
total J /ifj Xp distribution is more complex since it includes the feeddown from the ip' 
and Xc states. Then |^ 
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(13) 



Our results are shown in Fig. ^. We have chosen the octet absorption cross section 
such that the total J/'0 ot agrees in magnitude with the recent measurement by E866 
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Figure 4: The A dependence for color singlet and color octet absorption in the 
NRQCD model is shown. The results are calculated at 158 GeV (a), 450 GeV (b), 
and 920 GeV (c). The total J j^^ (solid), the direct J/-?/' (dashed), the ■?/'' (dot-dashed) 
and the Xc (dotted) A dependencies are given. 



at 800 GeV 



In this case, 0"°*^* 



5 mb and a^Jj^^ 



'j/i,N - ^ '^'^ ^^'^ ^j/^N - 2-5 mb gives a ^ 0.95 
at xj? ~ 0. The same octet cross section is then used for the octet component of 
the absorption for all charmonium states. If the octet state neutralizes its color, the 
resulting color singlet is absorbed according to Eq. (|10D- The same absorption cross 
sections are also used for the singlet component of charmonium production in section 
3.1. 

In Figs. and ^, the predicted difference between J/ip, ip' and Xc absorption was 
not large in the measurable region, particularly when the total J/ip A dependence 
was considered. Now, however, the Xc and total J/ip results are significantly different 
and if the NRQCD model provides the right description of charmonium production, 
the measured Xc A dependence should thus be quite different from that of the J/ip. 
Note also that the ip' a is slightly lower than the total J/ip a at x^ ~ 0, in accordance 
with the E866 results [g. 

The direct J/ip A dependence is rather similar to the octet results shown in Fig. 
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due to its large octet component. The ip' has a larger overall singlet component but 
the singlet influence on the A dependence is rather weak. The main difference between 
direct J/ip and tp' production at xj? > is the larger a of the tp' due to the singlet 
component. However, the dominant color singlet component of Xc production leads 
to an almost linear A dependence for < x^ < 0.5 at 158 GeV and —0.25 < xp < 0.5 
at 920 GeV. The range of xp at which a ~ 1 is broader at higher energies because 
the singlet gg contribution grows larger with energy. Given the similarities between 
the pure color singlet and color octet results at large negative xp, it is difficult to 
disentangle the relative contributions for the combination of the two in this region. 
However, the differences at large positive xp are due to the change in the relative 
octet/singlet contributions. At large Xp, the qq component is more important. This 
octet piece causes the drop in a of the Xc at large xp while having little effect on the 
J/ip and if)' . Finally, we note that the total J /iIj A dependence in this calculation is 
quite similar to the ip' dependence, as already indicated by previous measurements 

It is clear that if this model is correct, both NA60 and HERA-B should have no 
difficulty observing substantial differences in the J /ip and Xc A dependence since the 
values of a are clearly different even at positive xp. Other effects such as nuclear 
shadowing and energy loss would be similar for the two resonances so that differences 
in absorption mechanisms would not be washed out. 

4 Conclusions 

We have calculated the nuclear dependence of total and direct J /ip, ip' and Xc due to 
absorption alone. We have studied absorption of pure color singlets and color octets 
in the context of the color evaporation model and a combination of octet and singlet 
production in nonrelativistic QCD. When considering charmonium production in a 
pure color state, as in the color evaporation model, we find little difference in the 
charmonium A dependencies in regions accessible to past experiments, in agreement 
with the J /ip and ip' measurements to date. However, when the Xc is considered, its 
large color singlet component results in a substantially different A dependence in the 
nonrelativistic QCD description. This difference should be easily detected by the two 
experiments that plan to measure Xc production, NA60 and HERA-B. Their results 
should quickly answer the question posed by the title of this paper. 
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